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We propose a spherically symmetric and anisotropic model for strange quark stars within 
the framework of MIT Bag model. Though the model is found to comply with all the 
physical requirements of a realistic star satisfying a strange matter equation of state 
(EOS), the estimated values the Bag constant for different strange star candidates like 
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1. Introduction 

In 1964, Gell-ManrP and Zweig^ had independently suggested that hadr< 



ons are 



composed of even more fundamental particles called quarks, a proposition which got 
experimental support later on. In the theoretical front, the quark matter hypothesis 
put forward by WitterP^, had prompted many investigators to investigate an entirely 
new class of compact astrophysical objects composed of strange quark matter called 
strange stars (see Weber^ for a recent review). As quarks are not seen as free 
particles, the quark confinement mechanism have been dealt with great details in 
QCD. In the MIT bag modeP for strange stars, the quark confinement has been 
assumed to be caused by a universal pressure J5, called the bag constant. Farhi and 



JafeeMJ and Alcock et had shown that for a stable strange quark matter the value 
of the bag constant should be B ~ 55 — 75 MeV fm -3 . 

At the backdrop of such theoretical developments, analytic model building of 
strange quark stars naturally plays a crucial role to understand the gravitational 
behaviour of strange stars. In this vein, from the perspective of classical GR, we 
propose here an analytic solution capable of describing realistic strange stars. In an 
earlier work, some of U9 21 had proposed a relativistic model for strange stars where 
the geometric model of Krori and BarusP was used to describe the interior space- 
time of strange stars. In this paper, we adopt the Finch and SkeeP^ model to develop 
the interior space-time of a strange star. We assume a spherically symmetric star 
composed of strange quark matter whose equation of state (EOS) is governed by the 
MIT bag model. We also assume anisotropy in the matter composition implying that 
the radial pressure (p r ) is not equal to the tangential pressure (pt)- Since the density 
of a strange star may exceed the nuclear matter densi ty it is expected that pressure 
at the interior of such would be anisotropic, in general ^ Anisotropy in a highly 
dense compact stellar object may occur for various reasons (e.g, existence of a solid 
core, phase transition, presence of electromagnetic field etc.) whose contributions 
in the stellar modelling can be incorporated by assuming the star to be anisotropic 
in generate In our construction, the assumption of anisotropy provides an extra 
degree of freedom to deal with the system. 

The strange star model, in this work, has been developed by considering the 
MIT bag model EOS and a particular ansatz for the metric function g rr proposed 
by Finch and Skea^. The solution found here has been shown to comply with all 
the physical requirements of a realistic star. When applied to some proposed strange 
star candidates, it has been found that the corresponding values of the bag constant 
lie on the higher side as compared to its acceptable range put forward by Farhi and 
JafeeP and Alcock et aP 1 . In fact, our results are consistent with the experimental 
results from CERN-SPS and RHIC, indicating a wider range of values of the bag 
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2. Interior solution: 



We assume that the interior space-time of a strange star is described by a spherically 
symmetric metric of the form 



ds 2 = -e^dt 2 + + dr 2 + r 2 (d6 2 + sm 2 8d0 2 ), 



(1) 



where, the metric function v(r) is yet to be determined. No te t hat the ansatz for 
the metric function g rr in (JTJ) was proposed by Finch and SkegP^to develop a viable 
model for a relativistic compact star. The t = constant hyper-surface of the metric 
is paraboloidal in nature which exhibits a departure from spherical geometry and 
the constant R is a curvature parameter which governs the geometry of the back 
ground space-time. 

We assume that the energy-momentum tensor for the strange matter filling the 
interior of the star has the standard form T]y = diag(p, — p r , —pt, ~Pt), where, p is 
the energy-density; p r and pt are the radial and transverse pressures, respectively. 
Einstein's field equations for the line-element ([1]), accordingly, are obtained as (we 
set G = c = 1) 



S;r/ ' = iK 3+ ^) 



8irp r 
8np t 



1 + T7T 



1+^TT 



r 

R 2 

2 



r 
R 2 



R 2 ) 
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27 



R 2 



R?) 



-2 r 



V r 



(2) 
(3) 
(4) 



Eqs. (HJ-Q constitute a system of four unknowns (p, p ri pt, v). By suitably choosing 
any one of these unknown parameters, the system may be solved. Our objective here 
is to develop a model for strange stars and, therefore, we assume the simplest form 
of the strange matter EOS having the form 



Pr 



(5) 



where, B is the bag constant. Substituting Eq. ([5]) in Eq. ([3]) and integrating, we 
determine the unknown metric function v in the form 



1 , , 9 8-kB r 2 I 

v=-HR>+r 2 )-—- 2 [2R 2 



AttB 



(6) 



where is an integration constant. 

Using Eq. © in Eqs. Q, the density and the two pressures are then obtained 
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as 



8np - 
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3. Physical behaviour of the model 

We note that the physical behaviour of the model depends on the constants R and 
B. We need to put appropriate bounds on these parameters so that the model can 
describe a realistic strange star. To this end, based on various physical requirements, 
let us now analyze the behaviour of the physical parameters. 

Let us assume that b be the radius of the star. Then, from Eq. (j8]), the central 
and surface densities are respectively obtained as 



We also have 



Pa 



Pb 



8-kR 2 
1 

8t:R 2 



}P_ 

R 2 



dp 
dr 

g(r = 0)=0, 
dr 



R 2 



47ri? 4 (l + f^ 3 



<0, 



dr 2 



(r = 0) 



4?ri? 4 



< 0. 



(10) 

(11) 



(12) 

(13) 
(14) 



Obviously, the density is maximum at the centre and it decreases radially outward. 
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Similarly, from Eq. (JSJ), we have 

dp r r(5 



—2 — < 0, (15) 
ui Ltirnryi. -t- ■^J' 3 

%(r = 0) = 0, (16) 



d p r 5 
~dr T ~ ~ 12itR a 



< 0, (17) 



which show that the radial pressure also decreases from the centre towards the 
boundary. Thus, the energy density and the radial pressure are well behaved in 
the interior of the stellar configuration. Variations of the energy-density and two 
pressures have been shown in Fig. (p} and ([2]), respectively. 
The anisotropic parameter A(r) = - (p t — p r ) is obtained as 



A = — I — 



256tt 2 B 2 R 2 / r \6 / 5l2ir 2 B 2 R 2 64irB 32ir B \f r Y 
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Figure ([3]) shows the nature of the anisotropic stress at the stellar interior for a 
particular case. The condition that the anisotropic parameter A should vanish at 
the centre (r = 0) yields 

1 4 

— -B = 0, (19) 



8nR 2 3 

which can be used to calculate the Bag constant B. 



3.1. Matching conditions 

At the boundary of the star r — b, the interior metric should be matched to the 
Schwarzschild exterior metric. Continuity of the metric functions across the bound- 
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Fig. 1. Variation of the energy-density (p) at the interior of the star. 



p t 0.0006- 




~Pt PA 



Fig. 2. Variation of the radial(p r ) and transverse pressure(pt) at the interior of the star. 

ary surface yields 
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Fig. 3. Variation of the anisotropic stress A(r) = — (pt — p r ) at the stellar interior. 



From Eqs. (f20|). we obtain the compactness of the star as 



(22) 



The condition that the radial pressure must vanish at the boundary (p r (r 
b) = 0) determines the bag constant in the form 



D 



(23) 



3.2. TOV equation 

For an anisotropic fluid distribution, the generalized Tolman-Oppenheimer-VolkofF 
(TOV) equation gets the form 



dp r 1 ; / 2 , 

-T + o U (P + + ~ (Pr - Pt) = 0. 
dr I r 



We rewrite the above equation in the form 



M G (p + p r ) h=» dp r 2 

5 e 2 -j 1 — \Pt — Pr) = 0) 



(24) 



(25) 



where, Mc{r) is the effective gravitational mass inside a sphere of radius r and is 
given by 



M G (r) = -r^e^v'. 



(26) 
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Fig. 4. Three different forces acting on the fluid elements of a given configuration. 

The modified TOV equation describes the equilibrium condition for the strange star 
subject to gravitational (F g ), hydrostatic (Fh) and anisotropic (F a ) stresses within 
the stellar interior so that 

F g + F h + F a =0, (27) 
where the stress components are given by 

F g = -Br(p + p r ), (28) 

F a = - (pt - Pr) ■ (30) 
r 

In Fig. SI we have shown variations of F g , Fh and F a for a particular stellar config- 
uration. Obviously, it is possible to construct a static equilibrium configuration in 
the presence of anisotropic, gravitational and hydrostatic stresses. 

3.3. Energy conditions 

Imposition of the energy conditions, namely, the null (NEC), weak (WEC), strong 
(SEC) and dominant (DEC) energy conditions, put certain bounds on the model 
parameters. In our model, applying these energy conditions at the centre (r = 0), 
we get the following bounds: 

(i) NEC: po + po > => B < i.e. B < p , 

(ii) WEC: po + po > B < i.e. B < p , Po > 0, 

(iii) SEC: po + po > B < p , 3 Po + Po > =► B < f , 

(iv) DEC: po > \po\ =>■ B < f . 
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[ Strange Star - 4 U1S20 - 3ft R =7.117] 



Fig. 5. Variation of M e ff against radial parameter r. 

Values of the model parameters for different stellar configurations are in agreement 
with these bounds as shown in Table 2. 



3.4. Mass-Radius relation 

For a static spherically symmetric perfect fluid star in equilibrium, Buchdahff^ 
showed that the maximum allowed mass-radius ratio is given by 



271/ 



< 



(for a more 

generalized expression see Mak et aPj). In our model, the effective gravitational 
mass is obtained as 



M e ff = 4tt J pr 2 dr = 



b 2 
W 



b 2 
K 2 



(31) 



In Fig. [3 the effective mass for a given radius has been shown. The compactness of 
the star, accordingly, is given by 



M eff (b) 



1 



id 

R' 2 



(32) 



whose values at different radii have been shown in Fig. [6] We note that the constraint 
on the maximum allowed mass-radius ratio in our case is similar to the isotropic 
fluid sphere, i.e., A ' 



, < |. The corresponding surface redshift (Z s ) is obtained as 



[1 - (2u)Y 



1 



1 



1. 



(33) 



The maximum surface redshift for different strange star candidates in this model 
has been shown in Table 1. 
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[ Strange Star - 4 Ul 820 - 30, R - 7. 1 1 7] 




Fig. 6. Variation of — |££ against radial parameter r. 



Table 1 . Values of the maximum surface redshift for dif- 
ferent strange star candidates. 



Strange star candidate 


b (km) 


R (km) 


Zs(max) 


Her X-l 


7.7 


10.79 


0.2285 


SAX J 1808.4-3658(SS1) 


7.07 


5.787 


0.5787 


SAX J 1808.4-3658(SS2) 


6.35 


5.0282 


0.6108 


4U 1820-30 


10.0 


7.117 


0.7246 



3.5. Estimated bag values of some strange star candidates 

Based on the analytic model developed so far, to get an estimate of the range 
of various physical parameters, let us now consider some strange star candidates 
like 4U 1820-30, Her X-l and SAX J 1808.4-3658. Assuming the estimated mass 
and radius of these stars, we have calculated the values of the relevant physical 
parameters which have been compiled in Table 2. In Table 3, we have expressed 
the values of the bag constant, in particular, for respective stellar configurations. 
We observe that a wide range of values of the bag constant are possible unless we 
impose some (yet unknown) constraints on B. However, from the perspective of a 
mathematically self-consistent model, it appears that a wide range of values of the 
bag constant are possible which is consistent with the CERN-SPS and RHIC data. 

4. Conclusion 

We have obtained a new class of solutions for the interior of a compact stellar object 
like a strange star. To construct the model, we have used the phenomenological MIT 
bag model EOS for quark matter. The analytic solution obtained is non-singular 
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Table 2. Values of the model parameters for different Strange stars. (Data obtained for the star 4U 1820-30 have been 
utilized to plot figures.) 



Strange star candidate 


M (Mq) 


b (km) 


M 
b 


R (km) 


B (km~ 2 ) 


PO (km 2 ) 


Pb (km 2 ) 


PO (km 2 ) 


Her X-l 


0.88 


7.7 


0.168 


10.79 


0.00013 


0.00102 


0.00052 


0.000165 


SAX J 1808.4-3658(SS1) 


1.435 


7.07 


0.299 


5.787 


0.00021 


0.0035 


0.00085 


0.000901 


SAX J 1808.4-3658(SS2) 


1.323 


6.35 


0.308 


5.0282 


0.00026 


0.00472 


0.00107 


0.00121 


4U 1820-30 


2.25 


10.0 


0.332 


7.117 


0.0001 


0.0023 


0.00044 


0.00063 



Table 3. Physical values of energy density, pressure and Bag constant for different Strange stars. 



Strange star candidate 


Central density 
(gm cm -3 ) 


Surface density 
(gm cm -3 ) 


Central pressure 
(dyne cm -2 ) 


Bag constant 
(MeV fm- 3 ) 


Her X-l 
SAX J 1808.4-3658(SS1) 
SAX J 1808.4-3658(SS2) 
4U 1820-30 


1.38169 X 10 15 
4.808527413 x 10 15 
6.370476078 x 10 15 
3.179364432 x 10 15 


0.710012 X 10 15 
1.159209 x 10 15 
1.449092801 x 10 15 
0.5960016695 x 10 15 


4.145070271 x 10 35 
14.42558225 x 10 35 
19.11142825 x 10 35 
9.538093302 x 10 35 


99.7207866 
162.8102741 
203.5242699 
83.70809967 



and is anisotropic in nature. The estimated values of the bag constant B for different 
strange star candidates have been found to be on the higher side as compared to 
its estimated ran ge of 60 — 80 MeV fm -3 for a /3-equilibrium stable strange matter 
configuratiorPJ^J It is likely that a wide range of values of the bag constant are 
permissible which is in agreement with the recent CERN-SPS and RHIC data. 
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